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Abstract 

Using cosmological TV-body simulations of critical (SCDM) and open (il = 0.3, OCDM) cold dark 
matter models, we investigated the evolution of cluster galaxies. Based on our numerical simulation, we 
constructed merging history trees of the galaxies. By following their merging history, we could show that 
the major merger fractions of the galaxies in cluster-forming regions is roughly proportional to (1 -I- z)^'^ 
at low redshifts (z < 2), and has a steep peak at z ~ 2.5 and z ~ 3 in SCDM and OCDM, respectively. 
We also show that the cluster galaxies are affected by the tidal interaction after the clusters are formed. 
Because the formation redshift of the cluster in SCDM, Zform = 0.15, is much more recent than that of the 
cluster in OCDM, Zform = 1-6, the cluster galaxies in SCDM show more rapid evolution by tidal interactions 
from z = 0.5 than those in OCDM. 
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1. Introduction 

Rich clusters of galaxies are large laboratories for 
studying galaxy evolution (Dressier 1984), and their evo- 
lution can be followed with samples out to z ~ 1 (Rosati 
et al. 1998). It is well established that galaxy popu- 
lations vary with the density of neighboring galaxies in 
clusters of galaxies (Dressier 1980) and depend on the dis- 
tance from the clusters' centers (Whitmore et al. 1993). 
The increase in the fraction of blue, star-forming cluster 
galaxies with the redshift (Butcher, Oemler 1978, 1984a, 
1984b) has also been well established. It has been sug- 
gested that galaxy-galaxy and galaxy-cluster interactions 
play important roles in these effects; especially, a major 
merger of galaxies produces elliptical galaxies as merger 
remnants (Barnes 1989, 1996) and cumulative tidal inter- 
actions induce a morphological transformation of spiral 
galaxies to SO galaxies. 

Since it is generally believed that cold dark matter 
dominates the mass in the universe, we expect that the 
formation process of dark matter halos significantly af- 
fects the process of the formation and evolution of galax- 
ies. In this paper, we discuss our high-resolution cos- 
mological iV-body simulations which trace the motion of 
the dark matter particles and can resolve the galaxy-sized 
dark halos within a high-density environment. We also 
consider when the above mentioned interactions (i.e., ma- 
jor mergers and tidal interactions) act on the evolution of 
the cluster galaxies during the formation and evolution 
of clusters in a hierarchical clustering universe. 

We should consider the hydrodynamical processes of 
the baryonic components in order to follow the evolu- 



tion of galaxies. However, hydrodynamical simulations, 
e.g., smoothed particle hydrodynamics (SPH) simula- 
tions, need much more CPU time than coUisionless sim- 
ulations. Thus, it is difficult to obtain a wide dynamical 
range with such simulations. Here, we restrict ourselves 
to follow the evolution of dark matter halos and to use 
the galaxy tracing method described by Okamoto and 
Habe (1999, hereafter Paper I) to obtain merging history 
trees of galaxies. 

Recently, some authors have studied the evolution of 
the dark matter halos of the cluster galaxies in SCDM 
(Ghigna et al. 1998; Paper I). The epoch of the forma- 
tion of a cluster of galaxies is very sensitive to the value 
of the cosmological density parameter, Hq (Richstone et 
al. 1992). Since the evolution of clusters affects the evo- 
lution of galactic halos within the clusters (Paper I), it is 
interesting to compare the evolution of cluster galaxies in 
various cosmological models with different values of the 
density parameter. Here, we examine two cosmological 
models: one is the critical universe {Qq — 1), and the 
other is the open universe (JIq = 0.3). For both models 
we assume that the mass of the universe is dominated by 
cold dark matter (CDM). 

The plan of this paper is as follows: Techniques and 
parameters of the A^-body simulations and the method of 
creation of merging history trees of galaxies are described 
in section 2. Our results are presented in section 3 and 
discussed in section 4. 

2. Simulation 
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2.1. The Simulation Data Set 

Our simulations followed the evolution of a isolated 
spheres of a radius, i?sim, in both the standard CDM 
(SCDM) universe {Q = 1, h = Ho/WO kms-^Mpc"^ = 
0.5, as = 0.67) and the open CDM (OCDM) univ(-rsc 
(Oo = 0.3, h = 0.7, as = 1). The normalisations 
were chosen to approximately match the observed cluster 
abundance. Wc imposed the constraint of the 3a peak 
with an 8 Mpc Gaussian smoothed density field at the 
center of each simulation sphere in order to obtain a rich 
chist(!r (Hoffman, Ribak 1991). The simulations wore 
performed using a parallel tree-code, which was used in 
Paper I. To obtain a sufficient resolution to follow the 
evolution of the galaxy-sized halos with a relatively small 
number of particles, we used the multi-mass A^-body code 
(Navarro et al. 1997; Paper I). These initial conditions 
of our simulations were made as follows. 

First, only long- wavelength components were used for 
the realization of the initial perturbation in the simu- 
lation sphere using ^ 10''' particles; wc then performed 
a simulation with these low-resolution particles. After 
this procedure, we tagged the particles which are inside 
a sphere of radius 3 Mpc centered on the cluster center at 
z = 0. Next, we went back to the initial time stage, and 
then divided each tagged particles into 64 high-resolution 
particles according to the density perturbation that is 
realized by including additional shorter wavelength com- 
ponents. As a result, the total number of the particles 
became ~ 10^. 

We then calculated again the dark matter evolution us- 
ing high- and low-resolution particles from the new initial 
condition. Our analyses were operated only for the high- 
resolution particles. The mass of a high-resolution par- 
ticle was m ~ 5.5 X W^h~^M,^, and its softening length, 
e, was set to 5 kpc. 

The overall parameters and mass of the clusters in both 
simulations a,t z = are listed in table 1. 

2.2. Creation of Merging History Trees of Galaxies 

To create merging history trees of galaxies, we have to 
identify the galactic dark halos in the sea of dark matter. 
The identification of halos in such environments is a crit- 
ical step (Bertshinger, Gelb 1991; Summers et al. 1995). 
The most widely used halo-finding algorithm, called the 
friends-of- friends (e.g., Devis et al. 1985), is not accept- 
able, because it cannot separate substructures inside of 
large halos. Since the DENMAX algorithm (Bertshinger 
& Gelb 1991) shows good performance, we used its off- 
spring SKID (Governato et al. 1997). 

This algorithm groups particles by moving them along 
the density gradient to the local density maximum. The 
density field and the density gradient are defined ev- 
erywhere by smoothing each particle using the SPH-like 
method with the neighboring 64 particles. At a given red- 



shift, only particles with local densities greater than one- 
third of the virial density at that epoch are moved to the 
local density maximum. This threshold roughly corre- 
sponds to the local density at the virial radius. The final 
step of the process is to remove all particles that are not 
gravitationally bound to their parent halos. Here, halos 
which contain more particles than a threshold number, 
nth, are identified as galactic halos. Unless we explicitly 
state, nth = 30 is adopted. This is a large enough num- 
ber to inhibit the numerical evaporation of halos (Moore 
et al. 1996a). 

The method to create the merging history tree is simi- 
lar to that mentioned in Paper I; the details are as follows. 

We identify galactic halos with a 0.5 Gyr time interval, 
which is restricted by the disk space of our computer. 
Since this time interval is shorter than the dynamical 
time-scale of the clusters and the fading time-scale of 
evidence of starburst in galaxies, it is sufficiently short 
to construct the merger trees to investigate the effect of 
the merging of galaxies and the tidal interactions. 

The three most bound particles in each halo are tagged 
as tracers. We consider three cases to construct the 
merger tree of galaxies. 

First, if a halo at ij+i, where i is the number of a time 
stage, has two or more tracers that are contained in the 
same halo at ti, then the halo at is defined as a next 
halo of the halo at ti. In this case, the halo at <j is a 
progenitor of the halo at ti+i. 

Next, we consider the case that some halos at tj+i have 
one of three tracers of a halo at ti. The halo that has a 
tracer which is more bound in the halo at ti is chosen as 
the next halo of the halo at i^+i. 

Finally, if none of three tracers of a halo at ti are con- 
tained in any halos at the next time stage (ti), we refer 
the particle which is the most bound tracer of the halo at 
ti as a stripped tracer. We thus call both the halos and 
stripped tracers galaxies throughout this paper. 

We construct the merging history trees of galaxies in 
this way. A halo which has more than two progenitors at 
former time stage is referred as a merger. It often hap- 
pens that satellite galaxies pass through a central halo. 
Such cases should not be considered as merging. Hence, 
we check that the tracers of galaxies at ti-i, which are 
contained in a merger at ti are still in the same halo at 
ti+i- 

In order to estimate the mass of a stellar component 

of a galaxy, we assume that the mass of the stellar com- 
ponent is proportional to the sum of the masses of its all 
progenitors (hereafter we call this mass the summed-up- 
mass) . Except for the case in which a large fraction of the 
stellar component of the galaxies have been stripped dur- 
ing the halo stripping, this assumption may be valid. We 
estimate the summed-up-mass of a galaxy as the sum of 
the summed-up-masses of its all progenitors at the pre- 
vious time stage. For a newly forming halo which has 
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Table 1. Parameters of simulations. 
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Subscripts h and 1 indicate high-resolution and low-resolution particles, respectively. 




Fig. 1.. Mass evolution of the most massive virialized object 
in each model. The plus signs and diamonds represent the 
mass evolution of the SCDM cluster and the OCDM cluster, 
respectively. 

no progenitors at a previous time stage, the summed-up- 
mass of the halo is set to the mass of the halo. To consider 
mass increase due to accretion of dark matter to the halo 
after its first identification, we replace the summed-up- 
mass with the halo mass when the summed-up-mass is 
smaller than the halo mass. 

3. Results 

3.1. Evolution of the Cluster 

To define the size of clusters, we calculated the virial 
overdensity based on the spherical- collapse model. For 
the SCDM model we used 200 as the virial overdensity 
according to previous studies (e.g., Navarro et al. 1996). 
In OCDM, the virial overdensity is a function of redshifts. 
Therefore, we calculate it at each redshift and, it then 
became 6 ~ 400 at z = 0. In figure 1, we plot the mass 
evolution of the most massive virialized object in each 
model. It is well known that the formation epoch of the 
OCDM cluster is much earlier than that of the SCDM 
cluster. In our result, the formation redshift (the redshift 
when the half of the final mass has accreted) of the SCDM 
cluster is z ~ 0.15 and that of the OCDM cluster is 
z ~ 1.6. We show the x-y projection of a density map in 



a cube with sides of 2 ryir (?'vir is the radius of the sphere 
having the virial overdensity) centered on the cluster's 
center in each model at z = (figure 2). The gray scale 
represents logarithmic scaled density given by the SPH 
like method. We found that many galaxy-sized density 
peaks survive even in the central parts of the rich clusters. 

3.2. Merging of Galaxies 

Numerical simulations have shown that mergers with 
a mass ratio of 3:1 or less produce a remnant resem- 
bling an elliptical galaxy (Barnes 1996). Therefore, we 
define a halo which has more than two ancestors with 
this ratio at the former time stage as a "major merger." 
In figure 3, we show the major merger fraction of the 
large galaxies (lO^/i^^M® < Msum < 10^^/i"^Mg,) in 
the cluster- forming regions as a function of the redshift. 
The points in the figure can be fitted by a curve in 
Gottlober et al. (1999), a(l + 0)'^exp[7(l + z)], with 
a = 0.01,0.04, /3 = 4.5,4.6, and 7 = -1.4,-1.2 for 
SCDM and OCDM, respectively. When the clusters start 
forming, the merger fraction steeply decreases. One rea- 
son for this is that the high-velocity dispersion in clus- 
ters and groups inhibits the galaxies within these objects 
from merging with each other. Another reason is that the 
stripping of halos by tidal fields of such large objects pre- 
vents the merging of individual galactic halos (Funato et 
al. 1993; Bode et al. 1994). Since the cluster in OCDM 
forms much earlier than in SCDM (see figure 1), this de- 
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Fig. 3.. Major merger fractions of massive galaxies with 
Msum > 10^^h~^ Migi in cluster-forming regions. The plus signs 
and diamonds indicate the fraction in SCDM and in OCDM, 
respectively. 

clinc of the merger fraction appears at the higher redshift 
in OCDM than SCDM. 

Governato et al. (1999) showed that the major merger 
rate of the galaxy-sized dark halos in the field for z < 1 
is proportional to (1 + zf-"^ and (1 + zf-^ in SCDM and 
OCDM, respectively. To compare their result obtained 
in the field to our result obtained in the cluster-forming 
regions, we can say that the efficiency of the major merg- 
ing of galaxies in the cluster-forming regions more steeply 
increases toward high redshifts than in the field, espe- 
cially in OCDM. Since the density contrast of the cluster- 
forming regions in OCDM takes a larger value than in 
SCDM, the evolution of major merger rate in OCDM 
is significantly different between in the field and in the 
cluster-forming regions. 

Recently, van Dokkum et al. (1999) observed the rich 
cluster at 2; = 0.83, finding a high merger fraction in the 
cluster and rapid evolution of the fraction. The fraction 
in .2 < 1 is comparable to the result obtained here. 

3. 3. Tidal Stripping of Halos 

We evaluated the effects of the tidal stripping on the 
galactic halos in different cosmologies. For this purpose, 
we chose galactic halos at a redshift when the number 
of large halos with > lQ^^h~^M,^ is largest much 
before the cluster or group formation epochs. The tidal 
effects are probably negligible at such a redshift. Then, 
we examined whether the chosen galaxies loose their ha- 
los at lower redshifts. If they become Mh < 10^°/i~-^M® 
at lower redshifts, they must be tidally stripped. We 
then adopted nth = 19, 18 for SCDM and OCDM, respec- 
tively. Therefore, when a galaxy has become a stripped 



Fig. 4.. Stripped galaxy fractions of galaxies which have massive 
halos {Ml-, > 10"/i"^Mj) at the selected redshift in SCDM 
(left panel) and in OCDM (right panel). They are plotted in 
the 0.25h~^ Mpc bins from each cluster center at z = (plus 
signs) and at z = 0.5 (diamonds). 



tracer, it means that the galaxy does not have a halo 
with Mh > IQiO/i-^Mg,. 

When they become such small halos, dissipative effects, 
which were not included in our simulations, should be- 
come important. We call such galaxies stripped galaxies. 

In figure 4, we show the stripped galaxy fractions in 
the 0.25/i~^ Mpc radius bins from the cluster centers. 
The stripped galaxies show a similar distribution in both 
models at 2; = 0. Their evolutions, however, are very 
different between the models. In SCDM, there are few 
stripped galaxies, except for the central part at z = 0.5, 
and the fraction drastically increases near z = 0. On the 
other hand, we can see a strong correlation between the 
fraction and the radius in the OCDM cluster, even at 
z = 0.5; this fraction shows very weak evolution during 
z = 0.5 and 0. This is because the cluster in OCDM 
forms much earlier than in SCDM. 

Next, we compared the radii of the halos, which were 
determined as the radii at which their circular velocity 
profiles take minimum values (Ghigna et al. 1998; Pa- 
per I), to the tidal radii of the halos estimated at their 
pericentric positions, which we calculated by using the 
dark halo model of Navarro et al. (1997). The tidal radii 
of the halos, Test) were estimated by the following ap- 
proximation assuming an isothermal distribution of dark 
matter: 

^est — ^peri 77 ? V, / 

where Vmax is the maximum value of the circular veloc- 
ity of a galactic halo and Vc is the circular velocity of a 
cluster. In figure 5, we plot rest against rh for the outgo- 
ing halos that must have passed pericenter recently. We 
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Fig. 5.. Measured values of the radii of the member halos against 
the expected values, assuming that the halos have isothermal 
mass distributions and are tidally stripped at their pericentric 
positions. Only out-going large halos are plotted. The left 
panel and right panel represent the halos in the SCDM cluster 
and those in the OCDM cluster, respectively. 

plotted only the large halos {vc > 80 km s~^) to avoid 
any influence of the insufficient resolution. Moreover, we 
ignored the halos with Tperi < 300 kpc, because they have 
tidal tails due to impulsive collisions as they pass close 
to the cluster center (Ghigna et al. 1998; Paper I). In 
SCDM, most of the halos with rh > 100 kpc have larger 
radii than rest- On the other hand, the halos in OCDM 
are sufficiently truncated, and have comparable radii to 
Test, because few galaxies accrete to the cluster in OCDM 
at low redshifts (the cluster has formed at 2;form = 1-6). 
This means that the galaxies in the SCDM cluster tidally 
evolve, even at present, and such evolution has almost 
come to completion in the OCDM cluster. This result is 
consistent with the weak evolution of the stripped galaxy 
fraction of the cluster galaxies in OCDM. 

4. Discussion 

We have investigated the effects of the difference of the 
cosmological density parameter on the evolution of the 
cluster galaxies using the cosmological A^-body simula- 
tions in SCDM and OCDM. The cluster forms at ^! ~ 0.15 
and z ~ 1.6 in SCDM and OCDM, respectively. We have 
shown that the difference between the formation epochs 
of the clusters changes the evolution of the cluster galax- 
ies. 

The major merger fraction in the cluster-forming re- 
gions is roughly proportional to (1 + z)'^'^ in each cosmo- 
logical model at low redshifts (0 < z < 2) . The decline of 
this fraction appears at higher redshifts in OCDM than in 
SCDM. The reason is as follows. The efficiency of merg- 
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ing rapidly decreases as clusters form because of large 

internal velocities of the clusters and a reduction of the 
size of tidally truncated halos. Hence, the earlier forma- 
tion of the cluster in OCDM leads to an earlier decline of 
the major merger fraction. From this result, wc expect 
that in a lower density universe the elliptical galaxies in 
clusters would mainly form earlier. We will investigate 
this possibility in a forthcoming paper. 

The tidal interactions also have the possibility to 
change the morphology of the galaxies and to induce ac- 
tive star formation (Moore et al. 1996b, 1998). We find 
that, in the SCDM universe, the fraction of cluster galax- 
ies which have been stripped of their dark halos due to 
tidal interactions begins to increase from z ^ 0.5. Thus, 
if the morphological transformation from S to SO and the 
star burst due to the galaxy harassment (Moore et al. 
1996b, 1998) are caused by such tidal interactions, the 
morphology-density relation and the Butcher-Oemler ef- 
fect should evolve from z ~ 0.5 in SCDM. On the other 
hand, in OCDM, sice this fraction has already been sig- 
nificant at z = 0.5, we can observe these effects at higher 
redshifts than z = 0.5. For detailed analyses wc need a 
star-formation history in each galaxy, which will be con- 
sidered in future work. 

In this paper wc show how the cluster evolution affects 
on the evolution of cluster galaxies. The effects of differ- 
ent formation cpochcs of the clusters of galaxies between 
cosmological models on the color and morphological evo- 
lution of the cluster galaxies will be clarified by combin- 
ing the merging history trees obtained here with a simple 
model of gas cooling, starformation, and feedback used 
in semianalytical work (Kauffmann et al. 1993; Cole et 
al. 1994; Kauffmann et al. 1999). 

We wish to thank M. Fujimoto and M. Nagashima for 
useful discussions. Numerical computation in this work 
was carried on the HP Exemplar at the Yukawa Institute 
Computer Facility and on the SGI Origin 2000 at the Di- 
vision of Physics, Graduate School of Science, Hokkaido 
University. 
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